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Abstract 

Colloidal systems present exciting opportunities to study clusters. Unlike atomic clusters, 
which are frequently produced at extremely low density, colloidal clusters may interact with 
one another. Here we consider the effect of such interactions on the intra-cluster structure 
in simulations of colloidal cluster fluids. A sufficient increase in density leads to a higher 
population of clusters in the ground state. In other words, inter-cluster interactions perturb 
the intra-cluster behaviour, such that each cluster may no longer be considered as an isolated 
system. Conversely, for dilute, weakly interacting cluster fluids little dependence on colloid 
concentration is observed, and we thus argue that that it is reasonable to treat each cluster 
as an isolated system. 



1. Introduction 

Clusters are a distinct state of matter which exhibit different structural ordering and 
phase behaviour, relative to bulk materials Of particular relevance to, for example, 
many biological systems such as viruses, is their tendency to exhibit five- fold symmetry 
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SI] . Recently there has been a surge of interest in clusters 
110, 0, 0, G3, Q E3, Q, leading to the development 
These may in turn provide novel functionalised 



such as icosahedra and decahedra 
formed in colloidal systems 
of 'colloidal molecules' U 

materials fiol . E2I . ll4L Il5[. In addition to this technological relevance, visualisation of 
colloidal clusters allows direct access to their free energy landscape fl6| . 

Part of the attraction of studying colloidal dispersions is that, although in principle they 
are rather complex multicomponent systems, the spatial and dynamic asymmetry between 
the colloidal particles (10 nm-1 fim) and smaller molecular and ionic species has led to 
schemes where the smaller components are formally integrated out [13] • This leads to a one- 
component picture, where only the effective colloid-colloid interactions need be considered. 

Given that the structure of ground state clusters of simple liquid models is known, along 
with local energy minima [2], it seems natural to investigate the prevalence of such structures 
in colloidal systems, in particular those with depletion attractions such as colloid-polymer 
mixtures where, at fixed real temperature, an effective temperature may be interpreted as 
the inverse of the attraction strength between the colloids. These colloidal systems exhibit 
similar [18, 19, 2o| though not identical (l9, 21] structures to clusters of simple liquid models. 
In general, as the (effective) temperature is reduced, we expect more clusters in the ground 
state, unless kinetic frustration comes into play. 

Unlike atomic clusters, which are often considered in isolation, colloidal clusters may 
themselves form fluids 0, \^ . Colloidal cluster fluids are found in systems with competing 
interactions with short-ranged attractions and long-ranged repulsions 0, \^ . The attractions 
drive clustering while the repulsions prevent aggregation and phase separation, leading to 
a characteristic cluster size |22[. At sufficient concentration, these cluster fluids form gels 




arrest 



while for sufficiently strong repulsions, the cluster fluid can undergo dynamical 
or crystallisation [28|- Furthermore, the separation in length- and time-scales 
allows us to treat the system in a hierarchical manner. For example, dynamical arrest within 
clusters [29] and of a fluid comprised of clusters [26|, [27j can in general be decoupled. This 
hierarchy means that one may consider each cluster as an isolated system [30] . Alternatively, 
one may operate at the cluster-cluster level [31]. Here instead we consider the influence of 
the inter-cluster interactions on the mira-cluster behaviour. In other words, how is the free 
energy landscape of a given cluster perturbed by its neighbours? 

In a recent experimental study on a cluster fluid [30|, we observed that in fact a rather 
small number of clusters were found in the expected ground state, around 20%. However, 
after a careful mapping of experimental parameters to conventional spherically symmetric 
interactions, Brownian dynamics simulations showed that isolated clusters reached their 
ground state, unless geometric frustration led to ergodicity breaking for deep quenches [291 ]. 

Here we consider an experimentally relevant set of conditions [30] and use Brownian dy- 
namics simulation to investigate the effects of inter-cluster interactions on the intra-cluster 
behaviour. In particular we consider the population of clusters in the ground state. Our 
system is characterised by the colloid packing fraction 0, the strength of the attractive in- 
teraction /3sm and the strength of the repulsion /3ey, where f3 = ks is the Boltzmann 
constant and T is the temperature. We first determine the range of parameters that corre- 
spond to a cluster fluid, before investigating the effect that varying these parameters has on 
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the intra-cluster structure. 



2. Simulation details and model 

Although a great many improvements have since been made (32I . l33| . the theory of 
Asakura and Oosawa (AO) [34| is generally accepted to capture the essential behaviour of 
polymer-induced depletion attractions between colloids. This AO model ascribes an effective 
pair interaction between two colloidal hard spheres in a solution of ideal polymers. However, 
the hard core of the AO potential leads to difficulties with Brownian dynamics simulations. 
The Morse potential is a variable range spherically symmetric attractive interaction, and 
reads 



/3u M (r) = /3s M e^- r \e^-^ - 2), (1) 

where p is a range parameter, ^e M is the potential well depth and a is the diameter. 

Under the right conditions, this produces similar behaviour to the Asakura-Oosawa model 
(35| . Here we follow our previous work and set p = 33.06 and vary the potential well depth 
/Sem- This mimics the addition of polymer in the experiments and corresponds to a short- 
ranged attraction length of ~ 0.22 |29|. We used weakly charged colloids to suppress aggre- 
gation. Under some circumstances, one can treat the colloid-colloid electrostatic repulsion 
with a Yukawa form 

, N _ exp(— k(t — a)) / r\\ 

(3u Y (r) = (5e Y PV ) ^, 2 

r/a 

where k is the inverse Debye screening length. The contact potential is given by 

$ Ey T\~\ 7o\2~' (^) 

(1 + kg fT) 1 a 

where Z is the colloid charge and Ib is the Bjerrum length. We use the experimentally 
relevant value for the inverse Debye length kg = 0.5 and likewise consider two values for 
the contact potential of the Yukawa repulsion [3 By = 1,3 [30]. In the experimental system 
we seek to model, van der Waals attractions are largely absent. Colloids have a hard core, 
which is often added to equation El Here, however, the Morse potential provides a slightly 
softened core, so we consider the combined potential (3um + @Uy in the simulations. 

For the simulations we use a standard Brownian dynamics simulation scheme (36|. The 
scheme generates a discrete coordinate trajectory as follows 

D N 

Vi (t + 5t) = r i {t) + — - Fij(t)6t + fc?, (4) 

Kb-L - T~/. 

where St is the simulation time step, D is the diffusion constant and is the total number 
of particles. The colloids respond to the pairwise interactions and the solvent-induced 
thermal fluctuations Srf are treated as a Gaussian noise with the variance given by the 
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fluctuation-dissipation theorem. Treatment of the hydro dynamic interaction between col- 
loids is not included in the simulations. 

For each state point of fis Y = {1, 3}, (j) = {0.02, 0.05, 0.10} and fis M = {5, 6, 7, 8, 9, 10, 11, 12} 
which span the range from monomer fluids to aggregates/gels, we perform 8 statistically in- 
dependent simulations each of N = 1000 particles. Particles are initialized randomly subject 
to a non-overlap constraint in a cubic box. Periodic boundary conditions for the box are 
implemented. The Morse potential [eq. (Dp)] is truncated and shifted for r > 1.5a, where 
the Morse potential is typically less than 10 -6 . The electrostatic interactions are treated by 
adding a Yukawa repulsion term [eq. ([2D] for different values of (5ey as specified and this is 
also truncated and shifted for r > 10.175a, where the potential is of the order 10 -3 . We 
study the evolution of the system as the particles condense into small clusters under the 
influence of the Morse attractions, which are stabilized by the long range Yukawa charged 
repulsions. 

We define the Brownian time as the time taken for a colloid to diffuse its own radius: 

r (C7/2)2 ^ 



In the simulations, tb ~ 711 time units, while in the experiments tb ~ 9 s [30]. The time- 
step is St = 0.03 simulation time units and all runs are equilibrated for 1.5 x 10 7 steps and 
run for further 1 x 10 7 steps. The simulation runs therefore correspond to approximately 
1000 Brownian times or around 1 hour, a timescale certainly comparable to experimental 
work. 

We identify two particles as bonded if the separation of the particle centres is less than 
1.25a |29|]. Having identified the bond network, we use the Topological Cluster Classification 
(TCC) to determine the nature of the cluster (37j . This analysis identifies all the shortest 
path three, four and five membered rings in the bond network, and identifies clusters from 
these base units. We use the TCC to find clusters which are global energy minima of the 
Morse potential |2l| . We denote the number of particles in a cluster by m and follow Doye 
et.al. in terming the global energy minima clusters by the number of particles, the range 
of the potential and the point group symmetry of the cluster. The global energy minimum 
m = 3,4,5,6,7 clusters for the Morse potential with p = 33.06 are 3A triangle, 4A 
Td tetrahedron, 5A triangular bipyramid, 6A Oh octahedron, and 7A D$h pentagonal 
bipyramid respectively and the structures are depicted in Fig. [5l For m < 7 there is only 
one global minimum for all ranges of the Morse potential. For more details, see [371 ] . 



3. Results 

3.1. Phase diagram 

Here we are interested in cluster fluids, however at low density and/or low attraction, the 
system is dominated by unbound particles which we term a monomer fluid [Fig. Ha)] while 
at high density we find aggregation into much larger, often elongated [ist clusters and 
ultimately gelation [Fig. HTc)] pi. |23|. Between the aggregation/gel regime and the monomer 
fluid lies the cluster fluid seen in Fig. Hb). 
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Figure 1: Simulation snapshots of the three states seen (a) fiey = 3,0 = 0.02, /3em = 5, monomer fluid 
where more than 50% of the particles are identified as monomers/unbonded, (b) /fey = 3, = 0.02, f3sM = 9, 
cluster fluid in which we are interested, (c) /3sy = 1,0 = 0.10, fisM = 12, aggregate/gel where more than 
50% of particles are identified in m > 10 clusters. 




Figure 2: Cluster size distributions, m denotes number of particles in a cluster and lines are different values 
of fie M . (a) pe Y = 1,0 = 0.02, (b) f3e Y = 3, = 0.1. 
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Figure 3: Phase diagrams in the <\> — fisM plane. States are classified as monomer fluid, cluster fluid 
or aggregate/gel as defined in the text. The panels refer to different strengths of the Yukawa repulsion. 
Dashed lines are a guide to the eye. (a) /3ey = 1.0, (b) j3ey = 3.0. 



We begin our presentation of the results by determining those state points we consider 
to be cluster fluids. To do this, we analyse the cluster size distribution, as shown in Fig. El 
In the limit that ^Em ^ we expect a monomer fluid of isolated colloids as shown in Fig. 
[2](a) for weakly attracting systems f3sy — 1,0 = 0.02, ^e M ~ 4. We find a crossover to a 
system dominated by clusters at (3sm ~ 7. At higher interaction strengths, the number of 
monomers drops markedly, as the system moves towards larger clusters. As Fig. [2](a) shows, 
for moderate values of the attractive interaction f3sM, most clusters have a size m < 10. We 
take an arbitrary value of 50% to distinguish the three states in Fig. [TJ In other words, we 
consider a system with more than half the particles as monomers to be a monomer fluid, and 
with more than half the particles in m > 10 aggregates to be either aggregated or gelled. 
Increasing colloid packing fraction (j) = 0.1 and Yukawa interaction strength (5ey = 3 [Fig. 
[2](b)] leads to a somewhat different scenario. Here, we find clusters at all measured attraction 
strengths. Furthermore, for moderate interaction strengths, the distributions seem rather 
flat, with a high incidence of large clusters, or aggregates. In fact no monomers or dimers 
are present for /3sm > 8. Finding clusters at all interaction strengths is likely related to the 
increase in packing fraction to = 0.1 relative to Fig. E](a). Colloids have less chance to 
avoid one another. Recall the Yukawa repulsion is very long-ranged (kg = 0.5). The mean 
interparticle separation (60/7r) -1 / 3 at cj) = 0.1 is 1.74a, and the Yukawa repulsion varies 
only by 1.81/c^T between 1.74a and contact, so particles readily approach one another. 
Note that this crossover away from monomers at higher @£m is in marked contrast to a 
number of experiments 0, 0,0] which found an appreciable number of monomers even at 
high strengths of the attractive interaction. 

Based on the cluster size distributions, we now determine phase diagrams following our 
criteria for monomer fluids, cluster fluids and aggregates/gels. These are shown in Fig. [31 
The generic nature is similar to that found in previous simulation [23| and experimental 
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Figure 4: Pair correlation functions in cluster fluids, (a) weakly interacting cluster fluid, <p = 0.02, f3sy = 
1.0, Psm — 6.0. (b) strongly interacting cluster fluid, <\> — 0.1, /Sey = 3.0, fieM — 7.0 . 

0] work, of monomer fluids at low density and weak attractions, gels/aggregates at higher 
density /attract ion and a cluster fluid in between. Our definitions lead to a crossover rather 
than an abrupt transition for both the monomer - cluster fluid and cluster fluid - aggre- 
gate/gel. The differences between Figs. [3](a) and (b) are reasonably understood in terms 
of the increase in Yukawa repulsion, which inhibits clustering at the low colloid packing 
fraction (j) we consider. 

3.2. Inter- cluster interactions 

Pair correlation functions are shown in Fig. 01 Fig. 0](a) is a low-density (j) = 0.02 
fisy = 1.0 (3sm = 6.0 system, while for (b) (j) = 0.1 f3sy = 3.0 j3e m = 7.0. In both cases, the 
intra-cluster structure for r < 3a shows strong peaks. The pair correlation functions also 
allow us to comment on the strength of the cluster-cluster interactions. Fig. 0](b) shows a 
broad peak at r ~ 4a, conversely at longer ranges in Fig. 0](a), g(r) is almost flat. The peak 
in Fig. [H(b) indicates that these clusters exhibit strong correlations with their neighbours 
and are thus interacting with one another. Conversely, Fig. S](a) is a weakly interacting 
system. 

By thinking of the cluster fluid as a fluid of (polydisperse) Yukawa-like particles |3l| , we 
can gain some idea of the interaction between clusters. For the state point in Fig. [H(b), 
the mean cluster size is < m >^ 8.75. Typical cluster-cluster separations are given by the 
cluster-cluster peak i.e. rcc ~ 4a. If we consider clusters of average size at the typical 
separation, we have a cluster-cluster interaction (i.e. pair interaction) of f3Ucc —< m > 2 
(3Uy(rcc) ~ 14. Conversely, in the case of Fig. [H(a), there is a very weak peak around 
r « 6.25a, < m >~ 3.55, which leads to a typical cluster-cluster interaction strength 
of f3Uy(rcc) ~ 0.15. We can thus speak of weakly interacting [Fig. 01(a)] and strongly 
interacting [Fig. [H(b)] cluster fluids. We expect that a weakly interacting cluster fluid 
should show little change relative to isolated clusters, but that in a strongly interacting 
cluster fluid, the energy landscape of each cluster can be perturbed by its neighbours. 



7 




K K K 

Figure 5: Cluster populations of each cluster size 3 < m < 7 as a function of the well depth of the attractive 
interaction /3sm- Here we fix fiey = 3. We define cluster population by the ratio N c /N m where N m and 
N c are the total number of clusters of size m and the total number of clusters of type c (with m colloids) 
respectively. Symbols refer to colloid packing fraction, <j) = 0.02, circles, <j) = 0.05, triangles, (j) = 0.1, 
squares. The different plots corresponding to cluster structures (a) 3A Dsn triangles, (b) 4A tetrahedra, 
(c) 5 A Dsh triangular bipyramids, (d) 6Z C^v, (e) 6 A Oh octahedra, (f) 7A D$h pentagonal bipyramids. 
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Figure 6: Cluster populations for /3ey = 1. Symbols refer to colloid packing fraction, <\> = 0.02, circles, 
<j) = 0.05, triangles, <j) = 0.1, squares. The different plots corresponding to cluster structures (a) 3A D^h 
triangles, (b) 4A Td tetrahedra, (c) 5 A Dsn triangular bipyramids, (d) 6Z C2v> (e) 6A Oh octahedra, (f) 7A 
D$h pentagonal bipyramids. 
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3.3. Cluster populations 

In the cluster fluid, we treat each cluster individually, and measure its structure. In order 
to do this we identify the proportion of m-membered clusters which are in the ground state. 
We expect that, in the absence of geometric frustration (which is the case for m < 6), the 
cluster population will be dominated by the ground state at sufficient attraction. 

Our main results are given in Figs. [5] and El Here we see that, as expected, for 3 < m < 5, 
all clusters at all state points are largely found in the ground state for sufficient values of /3sm- 
The case of m = 6 [Figs. E](d,e) and[6](d,e)] is somewhat different, as two structures, the 6Z 
C2 V and the 6A Oh octahedron compete. Both have the same number of near-neighbours. 
Without Yukawa repulsions the difference in potential energy for our system is very small 
(of order 1 part in 1 x 10 8 ) [29], which slightly favours the 6A Oh octahedron. However 
symmetry and vibrational contributions to the free energy favour C<i v by a factor of around 
30 |16| as we indeed see. Here, the inclusion of Yukawa repulsions fisy > means that 6Z 
Civ is somewhat favoured energetically. 

In the case of 7-membered clusters [Fig. 0(f)] , the population of the global minimum 
7A D$h pentagonal bipyramid does not tend to unity. We showed this is due to ergodicity 
breaking at relatively deep quenches, such that geometric frustration can prevent access 
to the pentagonal bipyramid ground state [29]. For these relatively short simulation runs 
(r = 400t#), once the attractive interaction strength exceeds (5sm ^5.0 the bond lifetime 
exceeds the simulation run time. 

We now turn our attention to the role of increasing density for the fisy = 3.0 system 
(Fig. [5j). We see a general trend of higher density promoting access to the ground state. 
This is particularly true for those clusters in which there is no geometric frustration limiting 
access to the ground state, 3A Dsn triangles, 4A tetrahedra, 5A D 3 h triangular bipyramids 
and 6Z C<iv- The effect of colloid packing fraction appears strongest around ($Sm ~ 6. For 
example, we find an order of magnitude increase in the 6Z C^ v population at ($Sm = 6.0 upon 
increasing the packing fraction from (j) = 0.05 to = 0.1 [Fig. 113(e)]. For these clusters, the 
effect of raising the colloid concentration may thus be thought of as acting in a similar way to 
an increase in attraction Psm- The 6 A Oh octahedron population is around 1/30 that of the 
6Z Civ, which is consistent with simulations (29| and experiments (l6| on isolated clusters. 
Compared to isolated systems |29|, we see relatively few 7A D$h pentagonal bipyramids. 

In the case of weaker Yukawa repulsions (Fig. [6j), our arguments above concerning 
'strongly' and 'weakly' interacting cluster fluids would lead us to imagine that reducing the 
repulsions would lead to less inter-cluster interactions and less perturbation of the intra- 
cluster behaviour. This appears to be the case: compared to Fig. [5j the relatively weakly 
interacting clusters in Fig. [6] show less response to increasing the colloid packing fraction. 
In particular the more dilute (f) = 0.02 and 0.05 show little deviation from one another. 
This suggests that for these parameters, the cluster fluid approaches the dilute limit and 
behaves as a 'cluster gas'. Note that the experimental system considered in |30| was close 
to (j) = 0.02, fisy = 1.0 i.e. 'weakly interacting'. 
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4. Discussion 



We have considered the effect of cluster-cluster interactions on the the intra-cluster struc- 
ture in Brownian dynamics simulations of a cluster fluid. The overall behaviour is broadly 
similar to isolated clusters. That is to say, upon increasing the strength of the attraction 
clusters are able to reach their ground states, unless geometrically frustrated from doing 
so [29[] . For weak cluster-cluster interactions around 0.15/c^T between individual clusters, 
the intra-cluster behaviour depends little upon density. We therefore conclude that weakly 
interacting colloidal cluster fluids at low density (such as (j) = 0.02, fisy = 1) may be reason- 
ably treated as independent systems. Introducing a stronger coupling between the clusters 
such as increasing density leads to a higher population of clusters in the ground state for a 
given interaction strength. In other words, the energy landscape of each cluster is perturbed 
by its neighbours. 

This may be qualitatively understood in terms of the repulsive Yukawa interactions 
between clusters. The pair interaction energy between clusters typically increases as the 
clusters approach one another upon raising the colloid packing fraction. These inter-cluster 
repulsions might be expected to favour 'compact' clusters. One example of 'compact' is 
a cluster which minimises its radius of gyration. For spheres, these are 3A triangles, 
4A Td tetrahedra, 6A Oh octahedra and 7A pentagonal bipyramids (38|. Therefore, 
increasing cluster-cluster repulsions can lead to similar behaviour as increasing the strength 
of attraction. This may be equivalent to noting that clusters of hard spheres which minimize 
the radius of gyration exhibit the same structures as for the ground states for particles with 
attractive interactions 3A D^h triangles, 4A tetrahedra, 5A D^h triangular bipyramids, 
6A Oh octahedra and 7A D$h pentagonal bipyramids [4], One interesting approach would 
be to take a system which forms markedly different clusters, such as 'patchy particles' which 
can be tailored to form less compact structures in the ground state [10]. The competition 
between compact clusters favoured by density and the less compact ground state clusters 
could then be investigated. In this system, there in fact may be some competition between 
6Z Civ (favoured in isolation) and the 6 A Oh octahedron which is more compact. 

It is important to note that, although we simulate on experimentally relevant timescales, 
the cluster fluids we have studied here may not be truly at equilibrium. Colloidal cluster 
systems have been shown in simulation to form cluster crystals j28] or lanes (j^ j. but this 
suggests a more uniform cluster size distribution than we find here. Therefore, even though 
our cluster-cluster interactions can reach 14/cgT, the inherent polydispersity in the system 
prevents crystallisation. Systems with competing interactions often exhibit some form of 
structural ordering j40|]. It is reasonable to suppose that much longer simulation runs than 
we have been able to perform here might result in a shift in the cluster size distribution to 
one which is more monodisperse. Moreover, state points at higher concentration might be 



expected to develop into columnar phases or lamellae with a clear periodicity |40|, |4jJ]. We 



have not investigated the cluster dynamics in this system, but it is possible that there is a 
cluster glass region of the phase diagram. 

According to the definitions we have used here, there is no sharp transition between 
cluster fluids and aggregates/gels. Increasing density/attractions leads to aggregation and 
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gelation. However, locally in the gel we may expect the structure to resemble that of the 
clusters. This we indeed found in experiments on gelation without long-ranged repulsion 



42 



With competing interactions, as is the case here, a crossover from a cluster glassy state 
to a gel has been seen [27], where gelation is interpreted as a percolation phenomenon. That 
study found little change in local structure upon gelation. Thus both arrested spinodal type 
gels without long-ranged repulsions 14 21 [43 1 and those gels resulting from percolation in a 
system with competing interactions (y, |24|, |27|] may have a local structure dominated by a 



topology which follows that of isolated clusters. 



5. Conclusions 

The effect of cluster-cluster interactions on the intra-cluster structure of a model colloidal 
system has been studied. In the case of weak cluster-cluster interactions at low colloid 
concentration, the yield of structures is similar to if the clusters were isolated. This is the 
cluster gasElimit and interactions between clusters may be neglected. By increasing colloid 
density and/or increasing the strength of the Yukawa repulsion between colloids, cluster- 
cluster interactions become stronger and may no longer be neglected. Increasing either the 
density or Yukawa repulsion causes a higher yield of ground state structures, showing that 
energy landscape of a cluster is perturbed by the presence of neighbouring clusters. The 
cluster fluidEcannot be approximated as an isolated system due to the presence of cluster- 
cluster interactions which perturb the intra-cluster structure. 
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